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ABSTRACT 

We introduce Dirac processes, using Dirac delta functions, for short-rate-type pricing of financial 
derivatives. Dirac processes add spikes to the existing building blocks of diffusions and jumps. 
Dirac processes are Generalized Processes, which have not been used directly before because the 
dollar value of non-Real numbers is meaningless. However, short-rate pricing is based on integrals so 
Dirac processes are natural. This integration directly implies that jumps are redundant whilst Dirac 
processes expand expressivity of short-rate approaches. Practically, we demonstrate that Dirac 
processes enable high implied volatility for CDS swaptions that has been otherwise problematic in 
hazard rate setups. 

JEL classification: G12, C63. 
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Here we introduce Dirac processes, built from sequences of Dirac delta functions (Dirac, 1926 


Hoskins, 2009; Duistermaat and Kolk, 2010), and apply them to short-rate-type pricing of financial 


derivatives. Dirac processes expand the set of building blocks in mathematical finance beyond 
diffusions and jump processed built from Weiner and Poisson processes, to additionally include 
spikes based on the Dirac delta function. Dirac processes are a subset of Generalized Processes 


which were first introduced by Ito 

Ito, 

1954 

Gel’fand 

1955 

) some ten years after his introduction 

of Ito’s Lemma (Ito 

1944 

Doblin 

1940 

), but not previously used directly in mathematical finance. 


Generalized Processes have not previously been used directly in mathematical finance because the 
dollar value of values that are not in M, for example the value of a Dirac delta function at the origin, 
is meaningless. However, because short-rate pricing is based on integrals, Generalized Processes, 
and in particular the Dirac processes we introduce here, are natural. Dirac processes solve the 
general problem in short-rate pricing of how to introduce jumps in derived forward rates. This 
provides a required expressivity to short-rate-type modelling that puts it on a more equal footing 


with Forward-rate-type modelling which has used jumps for a long time ( 

Cont and Tankov 

2003 

Eberlein, Jacod, and Schonbucher 

2006 

Jiang and Yan 

2009 

Crosby 

2008 

Peng and Kou 

2008) 


For a particular example we show how to obtain high implied-volatility in options on credit default 


swaps (GDS swaptions) that is otherwise problematic in hazard-rate models (Jamshidian, 2004 


Brigo and Mercurio[ |2006t Kokholm and Nicolato, 2010[ Brigo and El-Bachir 2010; Roti, 2013[ 


Weekend 

2014; 

Stamm, Gallagher, and Lichters, 

2015 


Pre-crisis interest rate modelling was driven by the need to price hybrids and exotics leading 


to an emphasis on tenor Forward rate models such as SABR (Hagan, Kumar, Lesniewski, and 


Woodward 2002) and Levy-based models (jEberlein et ah 2006). In the credit space the emphasis 


was on structured credit (Brigo, Pallavicini, and Torresetti, 2010) and forward-based modelling 
hardly developed (Peng and Kou, 2008). Single-name credit modelling stayed with short-rate 


approaches for the most part (Brigo and Mercurio, 2006). Post-crisis, funding and credit costs, and 
regulatory capital costs, have driven the rapid development of valuation adjustments collectively 
called XVA (Kenyon and Stamm, 2012; Green, Kenyon, and Dennis, 2014; Green, 2015[ Stamm 


et ah, 2015), and XVA desks, as part of the need to price and manage counterparty netting sets 


together, mostly based on Monte Garlo techniques. This has brought short-rate and HJM modelling 
back into focus for interest rates. The Dirac processes we develop here fundamentally expand the 
expressivity of the short-rate and hazard rate approaches by delivering jumps in forward rates and 
enabling controllable smiles. They are equally applicable to commodities. 

Using a Dirac process for the hazard rate in GDS pricing means that there can be jumps in 
the survival probability of the reference entity from one day to the next. Given that the reference 
entity will generally only default from failing to make payments this provides an explicit model. For 
example, as in Figure it may have to redeem bonds on given dates, i.e. pay back the notionals, 
but have no comparable payments in the date range. Practically speaking the reference entity 


^Note that mean reverting processes are neither Levy processes nor Sato processes | Applebaum 

2004 

Sato 

2013 

but can be built from them (Schoutens and Carriboni 

2009 Kokholm and Nicolato 

2010 

I. 
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Figure 1. Example reference entity: known bond interest and principal repayments and possible 
refinancing repayments. Default may be unlikely outside of the specific redemption dates (or the 
previous weekends). 


may default on the prior Saturday to each payment, on the dates themselves, but have no other 
significant possibilities of default during the date range. Thus modelling hazard rates incorporating 
Dirac processes captures practically relevant dynamics. In addition, whilst CDS traders may focus 
on standardised CDS contracts that all, by definition, mature on only four dates per year (the IMM 
dates, Markit (2009)), credit pricing in general, e.g. for CVA or XVA (]Kenyon and Stamm 2012 


Stamm et ah, 2015; Green, 2015), requires arbitrary dates. 


Econometrically, jumps are observed in interest rates ( 

Das 

2002 Johannes! 2004 Piazzesi 

2005) 

and Forward-rate models with jumps are common ( 

Jiang and Yan 

2009). There is a wide literature 


on pricing interest rate options with jumps (Cont and Tankov, 2003) and general (Levy) approaches 
have been productive (Eberlein et ah, 2006) in the Eorward-rate space. However, without Dirac 
processes, jumps in forward rates or swaption prices cannot be created in the short-rate setting. 

Most stochastic processes used in derivative pricing are built from combinations of Brow¬ 
nian motions (continuous) and jump processes (discontinuous, i.e. Poisson processes). These 


can be characterized as good integrators which the Bichteler-Dellacherie theorem (Protter, 2010 


Bichteler, 2011) demonstrates are equivalent to semimartingales. Dirac processes, however, escape 


the Bichteler-Dellacherie theorem, and the definition of semimartingale, because whilst Dirac pro¬ 
cesses are adapted, they are not cadlag, nor do they take values exclusively in M. Instead, Dirac 
processes are ledc, (limites egale deux cotes), whilst mixed Dirac-Jump processes are, when at 
least one jump coincides with one Dirac delta function, lade (limites aux deux cotes). A related. 


but distinct, strand of research deals with functional extensions of Ito calculus (Cont and Fournie 


2010 ). 


Generalised functions have limitations, for example there is no single accepted method to deal 
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with multiplication of two generalised functions although methods exist (Colombeau 1992; Grosser 


1999). For our purposes we will always be using integrals of Dirac delta functions so this is not a 


significant limitation. We defer further theoretical positioning of Dirac processes, and their practical 
implications, until we have defined them. 

The first main contribution of this paper is the introduction of Dirac processes for short-rate 
pricing, expanding mathematical finance beyond diffusions and jumps to include spikes based on 
Dirac delta functions. The second main contribution is showing how Dirac processes can be used 
in derivative pricing including pricing CDS swaptions demonstrating high volatilities. A third 
contribution is the potential for detection of incomplete markets from single implied volatilities 
rather than an implied volatility smile. 

The paper is organized as follows: the first part introduces Dirac processes starting from their 
motivation in short-rate pricing, we then recall properties of the Dirac delta function; then we 
define the Dirac process and describe some of the (many) ways this can be used as a building block. 
The second part of the paper deals with interest rate derivative and single-name credit derivative 
pricing building up to pricing CDS swaptions with high implied-volatilities. Finally we discuss our 
results and conclude. 


I. Dirac Processes 

We first motivate our invention of Dirac processes starting from short-rate-type pricing, then 
recall features of the Dirac delta function, and then define the Dirac process and its relation to 
the broader category of Generalized Processes and stochastic integration. Given that the Dirac 
process is a building block in a similar way to Brownian motion or the Poisson process, there are 
an arbitrary number of ways to build with it. Equally, there is no need to use Dirac processes 
exclusively but they can be combined with diffusions and jumps. We introduce a basic set of 
constructions that is in no way exhaustive, but some of which are used in the following section on 
pricing derivatives with Dirac processes. 


A. Motivation 


Gonsider the price of a riskless zero coupon bond P{t,T), at t with maturity T, when the 
interest rate is modelled using the short-rate approach; 


P{t,T) = Ep 


g- 


( 1 ) 


r(t) is the short rate, i.e. the instantaneous spot interest rate. Qis the risk neutral measure which 
we assume exists, and E^ is the expectation with respect to the information available up to and 
including t (i.e. the appropriate filtration at t) and using Q. From Equation we see that only the 
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Figure 2. Processes (dashed) and their integrals (continuous line): (LHS) continuous process ; 
and (RHS) jump process. 


integral of the short rate is relevant, this is also the case for a defaultable bond P whose price is: 




e“ ft f{s)+X(s)ds 


( 2 ) 


where the hazard rate is given by A(t), the default time of the reference entity is r, and Ir>t is the 
indicator function for the survival of the reference entity. Again, only the integral of the hazard 
rate is relevant. In general observable and tradables only use integrals of short rates and hazard 
rate0 

Suppose we want to alter the volatility of these bond prices, we must alter the volatility of the 
short rate and the hazard rate. Consider Figure and the problem from an engineering point of 
view: what are the effects of introducing a jump process in the short rate on the integral of the 
short rate (or hazard rate) that determines the prices. We observe that there is zero significant 
qualitative effect from introducing a jump process: the integral of the short rate remains continuous. 
A non-differentiable point (without Generalized functions) is introduced, but if the short rate was 
based on a Brownian motion, which is nowhere differentiable (without Generalized Processes), then 
this is not a qualitative difference. Clearly, and (in retrospect) obviously, adding jumps to a short 
rate process is not useful from an engineering point of view to get qualitatively different volatility 
behaviour. 

If we add a unit shifted Dirac delta function to the continuous process in Figure we observe a 
jump in the integral as shown in Figure Thus we have achieved qualitatively different behaviour 
of the integral of the short rate. Thus a process based on Dirac delta functions is the appropriate 
engineering solution to constructing qualitatively different volatility in the integral. 

Ignoring the multi-curve nature of post-crisis interest rates for the moment (Kenyon, 2010| 


Mercurio 2010 Moreni and Pallavicini, 2014) consider a finite tenor forward interest rate from Ti 


^We sometimes use short rate for the approach, which encompasses any asset class, and sometimes for interest 
rates — this should be clear from the context. 


5 


























Value and Integral 



Time 


Figure 3. Continuous process with unit shifted Dirac delta function (dashed line, including vertical 
section at t = 1 extending to infinity but truncated in figure) and integral (continuous line) 


to T 2 defined as: 


1 


- !t^ r{s)ds _ 


( 3 ) 


T2-T1 

We see that Dirac delta functions in the short rate also create jumps in the forward rate which are 
otherwise impossible using only diffusions and jump processes. Now if we adopt a discount and 


tenor short-rate approach as in Kenyon (2010) (rather than the discount plus spread version in 


Kenyon and Stamm (2012)) the price of the Forward Rate Agreement is 


F{t,Ti,T2) = El 


_ r 
o Jt 


'^2 ^discount 


{s)ds_ 


1 


T2-Ti 


_ f T2 
3 Jti 


'^{s)ds 


- 1 


( 4 ) 


where “discount” indicates use of the discount-short rate and “tenor” indicates the use of the tenor 
discount rate, e.g. that for three month USD Xibor. Jumps in the Xibor rate dynamics can be 
present in the discount, the tenor, or both. 

A directly analogous setup for commodity forward rates using discount and convenience yields 


is possible in, say, oil forwards (Brigo, Chourdakis, and Bakkar, 2008) to introduce jumps there. 


Let us now turn to default timing. Consider Figure which shows an example of a known 
repayment schedule for bond principal and interest. In the short term, say up to six months in 
this example, there are many known repayments, but fewer thereafter. In general this pattern does 
not mean that the reference entity becomes debt-free but simply that the debt refinancing has not 
yet been arranged. Thus we have a mixture of deterministic times with stochastic event outcomes 
(default or not) and a region of stochastic event times and outcomes. Many reference entities share 
this pattern of known and uncertain events. The key point is that default is very unlikely outside 
of the repayment dates (and the previous weekends for non-sovereign borrowers) for the first six 
months. Most defaults of significantly-sized entities are liquidity defaults (i.e. failure to refinance) 
not insolvency defaults (declaration by auditors). This mix of event types can be modelled using a 
Dirac process with mixed deterministic and stochastic event arrivals. 
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Approach to a Dirac Delta Function Approach to a Heaviside Step Function 




Time Time 

Figure 4. Sequence of functions approaching a Dirac delta function (LHS) and its derivative the 
Heaviside function (RHS). 


B. The Dirac delta function 


We provide the basic details of Dirac delta functions here for convenience based on (Arfken 


2012; Hoskins, 2009; Duistermaat and Kolk, 2010). The Dirac delta function was first introduced 


in its modern form for quantum mechanics (Dirac, 1926) and is widely used in physics (partial 


differential equations, Green’s function, Fourier analysis, Arfken (2012)). A Dirac delta function 


is an example of a Generalized function (Hoskins, 2009), also known as distributions (Duistermaat 


and Kolk, 2010). A rigorous mathematical treatment based on distribution^ was built in the 1950s 
and extended to algebras (i.e. treating multiplication of generalized functions by each other) in the 
late 1980s ( Colombeau[ 1992; Grosser, 1999). 

Theoretically, generalized functions can be viewed as the smallest closure of the space of usual 
functions (i.e. those with values in M) under differentiation Duistermaat and Kolk (2010). They are 


not ordinary functions in that their properties can generally only be understood in the context of 
their effects on other functions, termed test functions, or upon integration to create usual functions. 
A Dirac delta function can be considered as a limit of a sequence of functions: 


fk{t -a) = 


1/k a < t < a + k 
0 otherwise 


The particular sequence of functions is not unique, many other sequences are possible (e.g. Figure 
, but the properties of the limit are unique in an appropriate sense (Hoskins 2009; Duistermaat 
[and Kolk 2010), for example with respect to their behaviour on a universe of test functions. 

Now the integral is constant, i.e.: 


ra-\-k 


h= j fk{t - a)dt = 

Jo 


1 

k^ 


dt = I 


( 5 ) 


^This re-uses the word from probability with some overlaps. 
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Thus the Dirac delta function can be defined as: 


6{t - a) = lim fk{t - a) (6) 

fc —^0 

Obviously this is not a usual function. Considered as a usual function Equation implies that the 
function is zero everywhere except at a where it is infinity, which is not an element of the Real 
numbers. 

Apart from the integral property in Equation the other key property we will make use of is 
its sifting property (with no “h”) on any ordinary function g{t): 

g{t)6{t - a)dt = g{a) (7) 

The usual distributive property of differentiation holds, i.e. 

(g(t)6(t - a))' = g{t)'5{t - a) + g{t)5{t - a)' 



where 

g{t)6{t - a)'dt = -g{a)' 

We will not consider multiplication of generalized functions as this is not needed for our applications. 

The Dirac delta function can also be considered as the differential of the Heaviside step function, 
i.e. formally 

^u{t - a) = S{t - a) = Sa{t) (8) 

where u{t) is the Heaviside step function, u{t) = 0,V t < 0; 1,V t > 0. 

Let us assue that g{t) is some well-behaved test function, where well-behaved means that it 
takes values in M and is bounded. Using the Riemann-Stieltjes integral we see that integrating this 
test function multiplied by a Dirac delta function is equivalent to integrating the test function with 
respect to a Heaviside step function, i.e. 

/ OO POO 

g{t)5{t — a)dt = / g(t)du{t — a) = g{a) 

-OO j —OO 



This provides another connection between Dirac processes and jump processes, and an 
definition of the Dirac delta function. That is, since integrals of Dirac delta functions 
behaved they have a representation in terms of ordinary functions, i.e. as step functions, 
see that the derivative of a Poisson process is a Dirac process and vice versa. 

A Dirac delta function can also be regarded as a functional, rather than a Generalized 


for its effects on test functions within integrals Hoskins (2009); Lighthill (2003). 


alternate 
are well- 
We shall 

function. 









C. Dirac processes 

We start by defining a Dirac process and then build upon it. Just as a Brownian motion alone 
is not particularly useful, similarly a Dirac process is useful for what can be done with it. 

DEFINITION 1 (Dirac process): A Dirac process D{t) is a sequence of Dirac delta functions (5(*) 
separated (time shifted) by exponentially distributed waiting times ej, i.e. 


OO 

D{t) ■.= Y,5{t-E,) 

i=0 

j=i 

Ei = ^ Cj 
3=0 

Cj ^i.i.d. Exponential random variable, mean v 

PROPOSITION 1 (Dirac process properties): A Dirac process is Markov, memoryless, and non- 
anticipative. 


Proof. Proof of Dirac process properties: Since the Dirac delta function exists only at a single time 
point, and the time between Dirac delta functions is exponentially distributed and the exponential 
distribution is memoryless all three properties follow. □ 


The proposition above relies on the fact that the Dirac delta function is the limit of the function 


sequence (as in Equation LB or similar) rather than any element of that limit. If this were not the 
case then there would be a potential issue with being able to predict some time in advance when 
the next spike would occur. However, this is not the case. 

The rate of arrival of spikes in a Dirac process is constant with rate u. 


DEFINITION 2 (Non-homogeneous (or inhomogeneous) Dirac process): A Non-homogeneous Dirac 
process has a rate parameter iy{t) that is a deterministic function of time. 

DEFINITION 3 (Compound Dirac process): A Compound Dirac process Dx{t) is a Dirac process 
where each spike is scaled by a non-anticipative stochastic process x{t, i) 


OO 

Dx{t) := ^ x{i, f) 6{t - Ei) 
i=0 

Note that the scaling process, x{i,t), above does not need to be Markov and does know which 
spike it is scaling (via i). The effect of the scaling process is to change the integral of the process 
as in Equation!^ We have also not ruled out that the scaling is affected by the time at which the 
spike occurs (t). 

DEFINITION 4 (Deterministic Dirac process): A Deterministic Dirac process is a sequence of 
time-shifted Dirac delta functions where the time-shifts are known in advance. 
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DEFINITION 5 (Compound Deterministic Dirac process): A Compound Deterministic Dirac pro¬ 
cess is Compound Dirac process where the time-shifts are known in advance. 


A Compound Deterministic Dirac process models a known sequence of events where there is 
uncertainty as to the magnitude of the events, but not when they should occur. In interest rate 
modelling, for example, the meeting dates of the Bank of England’s Monetary Policy Committee 
(MPC) are knowrj^ as are the meeting dates for the Federal Open Market Committee (FOMC)[^ 
at least a year in advance. Jumps have been associated with these meeting dates ( jPiazzesi] |2005[ ). 
In credit as in Figure repayment dates may be known but not the precise probability of default 
on each repayment date. In the longer term, for credit, there are repayments from new refinancing 
so both a Compound Deterministic Dirac process and a Compound Dirac process can be required. 
Equally the exact dates of MPC and FOMC committee meetings are not known more than a year 
in advance. 

Now that we have defined the Dirac process and some variants we will place it into context 


with respect to stochastic integration from the Protter (2010) and Bichteler (2011) perspective in 


general, and with respect to Poisson processes in particular with the two following theorems. We 
will then demonstrate with the third theorem just how flexible the Generalized process framework 
is. The first two theorems are new, whereas the third theorem is known but all three are useful for 
positioning Dirac processes. 

THEOREM 1 (Not good Bichteler-Dellacherie integrator): A Dirac process is not a good integrator 
in the Bichteler-Dellacherie (Protter. 2010: Bichteler. \2011) sense. 


Proof. Proof of Not good Bichteler-Dellacherie integrator: Obvious because it is not a semi¬ 
martingale because it is not a usual function, i.e. taking values only in M. □ 


THEOREM 2 (Integrator-integration equivalence): Integrating with respect to a Poisson process is 
identical to integrating a Dirac process with respect to time and vice versa. 

Proof. Proof of Integrator-integration equivalence: Direct consequence of Equation □ 

The equivalence theorem above shows that we have the choice, formally, of integrating with 
respect to a jump process or integrating a Dirac process w.r.t. time in terms of modelling. Generally 
we chose to integrate Dirac processes as the most direct and natural approach. This approach 
also avoids the potential confusion of specifying whether a jump process is the integrand (the 
function being integrated) or the integrator (the function that the integration is with respect to). 
This potential confusion is particularly acute in credit if we are comparing with survival process 
approaches where there are jumps in the survival process (which is the integral of the hazard rate 
process). 

^http://www.bankofengland.co.uk/publications/Pages/news/2014/119.aspx 
“http://www.federalreserve.gov/monetarypolicy/fomccalendars.htm 


10 



















Series approaching a Compound Dirac Process 



Integrated Compound Dirac Process 
aka Compound Poisson Process 



Figure 5. Approach to a Compound Dirac process and its integral the Compound Poisson process. 

THEOREM 3 (Integration w.r.t. a Dirac process): Integration with respect to a Dirac process 
defined for each Dirac delta function as 


g{t)d5{t — a) = / g{t)5'{t — a)dt = —g'{a) 


(9) 


is well defined. 

Proof. Integration w.r.t. a Dirac process: Elementary consequence of Dirac delta function proper¬ 
ties and closure of Generalized functions under differentiation. □ 

The three theorems above serve to define integration of Dirac processes and integration by 
Dirac processes with respect to standard stochastic integration theory (Protter, 2010| Bichteler 


2011). Although a Dirac process is not a good integrator it is possible to integrate with respect to 


it. We will only do so formally for convenience, but the machinery is well-defined via Generalized 
functions. Integration of, rather than by, Dirac processes is also illustrated by the relationship 
between a Gompound Dirac process and its integral the Gompound Poisson process as in Figure 
Qualitatively, integration of (not by) the respective processes moves rightwards below. 

Dirac process Jump process Gontinuous process 

It is possible to differentiate Brownian motion in the space of Generalized processes, but this is has 
not proved generally useful in mathematical finance although Generalized Processes are mentioned 


in the context of white noise Seydel (2012). Ito’s first approach (Ito, 1944) via Ito integrals is 


the dominant method in derivative pricing. Ito later invented Generalized processes (Ito, 1954 


Gel’fand, 1955) but these have not previously found general application in mathematical finance. 
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Figure 6. USD CDS market size from quotes from a major data provider pre-selected by the 
provider for at least minimal liquidity. 


II. Pricing with Dirac Processes 

In the previous section we defined the Dirac process and a set of variants, now we turn to 
pricing financial derivatives using these Dirac processes. We will start by considering deterministic 
Dirac (spike) arrival rates which are sufficient for linear products and then move to stochastic 
Dirac processes for non-linear products. We will price an option on a zero coupon bond and then 
a CDS swaption in a hazard rate setup. We pick CDS swaptions because they have previously 
been problematic to price in hazard rate setups even with jumps, i.e. it has been difficult to obtain 


sufficiently high implied volatilities (Jamshidian, 2004 Brigo and Mercurio 2006; Kokholm and 


Nicolato 2010; Brigo and El-Bachir, 2010; Roti 2013; Weekend, 2014; Stamm et ah, 2015). We 


stress that the models presented are the first examples of the use of Dirac processes. Just as there 
are an arbitrary number of models that can be built from diffusions and jump processes, the same 
is true of Dirac processes alone and in combination with diffusions and jump processes. 


A. CDS and CDS swaptions 

The USD CDS market coverage grew rapidly from around 200 reference entities in 2002 to 
around 1600 in 2008 and remained at roughly that size thereafter, see Figure Coverage is 
relatively small compared to the number of counterparties of any major bank. We interpret CDS 


spreads below in terms of protection and premium legs, as is standard (Brigo and Mercurio 2006) 


However, this ignores the fact that credit protection provides regulatory capital relief (BCBS-189 


2011), which may be significant in interpreting CDS spreads (Kenyon and Green 2013). We leave 


this extension dealing with incorporating capital pricing along the lines of (Green et ah, 2014) for 
future work. 

The single-name CDS swaption market is almost exclusively OTC (over-the-counter, i.e. be¬ 
spoke). Historical CDS spread volatility can be very high and level dependent, see Figure]^ The 
data in the figure are for reference entities with (almost) complete data series for 5Y CDS quotes 
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5Y CDS 2007-2011 



Figure 7. Density plot of historical 5Y CDS relative volatilities and spreads mid-2007 to end-2011 
using a one-year window for both. Plot combines data from 122 reference entities with at least 
99% data coverage of the period. CDS are senior (SNRFOR), USD-denominated, and with XR 
document clause (most common type). Volatility is in percent and Level is in bps for the mean 
CDS spread (each taken from overlapping one-year windows). This is a downward-biased sample, 
at least for CDS level, because of the requirement for (almost) complete data across the crisis. 
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from mid-2007 to end-2011 so there is a significant downward bias in the data as quotes dried up 
on many financial (and other) reference entities during the crisis. We do observe that it is certainly 
not true that lower CDS spread always implies lower CDS spread volatility, at least historically. 
Individual reference entities display such a wide variety of patterns that the summary across all 
reference entities is given here. Market-implied volatilities may be much higher than historical, 


especially considering CDS index volatilities, e.g. CDX Stamm et al. (2015) which have crossed 


100% several times. Note that here we focus on single-name credit and leave multi-name credit for 
future research. 

B. Short Rate and Hazard Rate Pricing 

Here we give the model-independent formulae that we will expand upon in later sections. For 
credit derivatives the short rate approach is usually called the hazard rate approach. 

The price of a riskless (i.e. non-defaultable) zero coupon bond is 


P(t,T) :=Et 


g- h riu)du 


where r{u) is the short rate. A defaultable zero coupon bond is 


P{t,T) :=Et 


g- r{u)+X{u)du 


where X{u) is the hazard rate. A semi-defaultable zero coupon bond is 


P{t,Ti,T2) :=Et 




The semi-defaultable case is where the default risk and the discounting have different maturities. 
This is not a (currently) traded instrument but is used below in CDS pricing and is the general 
case. 

A CDS contract provides protection against default of the reference entity in exchange for a 
series of periodic payments called coupons. Thus there is a protection leg paying 1-recovery on 
reference entity default, and a premium leg where coupons are paid provided that the reference 
entity has not yet defaulted. We will neglect accruals for simplicity, they can be added simply. 
Since the Big Ban^CDS spreads are quantized, i.e. they take only a certain set of standard values 
and any difference is made up using an up-front fee. Many other features are also standardized but 
these details are not material here. 

Default is not, practically, a continuous process but rather only happens with a minimum 
temporal resolution of about a day, so we use a discrete expression for the protection leg rather 


than the integral form that is more commonly seen (Brigo and Mercurio, 2006). This is more 


accurate and can be expressed in terms of the primitives we have already, i.e. defaultable bonds 
and semi-defaultable bonds. 


^CDS Market, not the astronomical one. 
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The formula below is model independent. 

THEOREM 4 (CDS Price): The price of a CDS for protection [Ta,Tf,] with premium R is 

CDSafi{R) = ProtectionLeg^ — PremiumLegg^i,{R) (10) 

d b 

= LGD Y, (R(0,T,)-P(0,T„ri_i))-R aiP{0,Ti) (11) 

2 =c+l i=a-\-l 

where Oi is the year fraction, and [Tc ,..., T^] gives a daily partition of [Tq, Tf\, and LGD is the loss 
given default. 

Proof. Proof: The Premium leg is the discounted coupon payments conditional on the reference 
entity not having defaulted. This is a sum of defaultable bonds scaled with the premium payments. 
The Protection leg is the sum of default risks over each day discounted back to time zero. □ 

C. Pricing with Deterministic Intensity and Severity 

This section gives prices for zero coupon bonds, riskless, defaultable, and semi-defaultable with 
deterministic intensity and severity. In the next Section we cover the stochastic case and option 
pricing. 

C.l. Zero Coupon Bonds 

Consider a general deterministic Dirac case where we use a differential shorthand for hazard 
rate A(f) given by the s-scaled inhomogeneous Dirac process 

dX{u) = s dD^^u){u) 

where s is the (deterministic) severity of each spike. This equation states that the hazard rate is 
zero except where it is a s-severity Dirac delta function, and the process is inhomogeneous with 
(rate) intensity C{u). It is clear that X{u) is independent of r{u) assuming that the interest rate 
short-rat^ has no jumps and no Dirac process component. 

THEOREM 5 (Bonds: Dirac Hazard with Continuous Short Rate): The price at time t of a de¬ 
faultable zero eoupon bond maturing at time T with unit notional, assuming 

• the short rate is continuous; 

• the hazard rate is given by dX{u) = s dDQ(^W}{'^) where s, ((u) are deterministic; 
is given by: 

P(t, T) = P(t, T) exp ^(e“^ — l) J f{u)du 

^We will use short-rate for interest rate short-rate hereafter, and hazard rate for default rate short rate. 


15 



and the price of a semi-defaultable zero coupon bond is given by: 

P(t,ri,T 2 ) = P(t,ri)exp - 1) Ciu)du 

Proof. Proof: We will prove the semi-defaultable case as it is the more general. Since there are no 
jumps in the riskless bond price the semi-defaultable bond price is independent of it, by indepen¬ 
dence of Dirac and Brownian processes. Thus: 


g- ft ^ Mu)du 


p(t,ri,T2) = p(t,ri)Ei 

= P{t, T 2 ) exp ^(e“^ ~ 


( 12 ) 

(13) 

(14) 


where ({u)du) is a random variable with a Poisson distribution with scale C{u)du. □ 


T2 


Implicit in these equations is the fact that the event arrival and severity can be traded off against 


each other. Making the approximation (for now and for clarity) taken from Brigo and Mercurio 
(2006) that 


hazard rate 
A 


CDS spread 
Loss given default 
CDS 
LCD 


(15) 

(16) 


then equating the two expressions for a defaultable zero coupon bond we have 


P(0,t)e-^* = P(0,t)e(" 


(17) 


which implies that the tradeoffs between arrival rate u and severity s, are 

A 


Figure shows the tradeoffs for two CDS levels. These tradeoffs will be significant for multi-name 
credit to change the default-time correlation but since this is out of scope for this paper we do not 
go further here. 

We now consider the case where the short rate is also a Dirac process. After this we will consider 
the mixed case where the short rate is a mix of continuous and Dirac processes. 

THEOREM 6 (Bonds: Dirac Hazard with Dirac Short Rate): The price at time t of a defaultable 
zero coupon bond maturing at time T with unit notional, assuming the short rate and hazard rates 
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Severity 


Figure 8. Tradeoffs possible between event arrival rate and severity for two CDS spread levels 
(200bps and 400bps) using Equation 
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see text for details. 


are: 


dr{u) = („)(«) 

d\{u) = S 2 d%(„)(M) + sodD(-^^^^{u) 


is given by: 


P(t,T)) =exp ^2 (e - l) ^ Co(«)drt+(e - l) ^ Ci{u)du 
+ - l) ^ C,2{u)dv^ 

and the priee of a semi-defaultable zero coupon bond is given by: 

P(t,T)) =exp - l) ^ Co{u)du + {e~^^ - l) Ci{u)du 

+ (e"*° “ ^ Co{u)du + (e "®2 - l) ^ C 2 {u)du^ 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


Proof. Proof: Almost surely there are no common spikes between different Dirac processes Cii^)- 
Now: 


P{t,T))=Et 


^-{soNiff Co(u)du)+Si_N(ff' ^i{u)du)+soN{ff' Coiu)du)+S2N{ff' (;2{u)du)) 


(25) 


each Poisson process N{*) is independent, so their exponentials are independent because functions 
of independent random variables remain independent so the result follows. The semi-defaultable 
bond is a direct extension. □ 


THEOREM 7 (Bonds: Dirac Hazard with Mixed Short Rate): The price at time t of a defaultable 
zero coupon bond maturing at time T with unit notional, assuming the short rate and hazard rates 
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are as in Theorem except that the riskless bond price is: 


P{t,T) = Pc{t,T)Pn{t,T) 


where Pc{t, T) is the part of the price given by a continuous short rate and Poit, T) the part of the 
price given by a Dirac process, is given by: 


P(t,r)) =Pc'(i,r)exp ^2 (e - l) ^ Coiu)du + {e - l) Ci{u)d' 

+ (e"®2 - 1) ^ C2iu)du^ 


and the price of a semi-defaultable zero coupon bond is given by: 

rP 


P{t,T)) =Pc{t,Ti)exp (^{e - 1) Co{u)du + {e - l) Ci{u)du 


+ (e - l) f Co{u)du + (e - l) f C. 2 {u)du 

Jt Jt 


(26) 

(27) 


(28) 

(29) 


Proof. Proof: There is no correlation between a continuous process and a Dirac process, so the 
result follows given Theorem]^ □ 

Calibration to yield and survival curves are direct. We observe, however, that we can trade 
off the severity against the intensity for the Dirac process but this has no effect on the pricing 
provided that there is no tradeoff between different Dirac processes. This could affect multi-name 
credit pricing but that is out of scope for this paper. 


C. 2. Interest Rate Swaps 

We confine our attention to interest rate swap (IRS) pricing using a single-curve which we 
take as riskless. Typical examples are SONIA, EONIA and Fed Funds curves for GBP, EUR, and 
USD respectively. Multi-curve pricing is outside the scope of this paper as the focus is on credit 
applications. 

Given riskless discount bond prices IRS prices are standard based on Equations and making 
the usual assumption of fully collateralized trades with cash collateral in the same currency. We 
neglect initial margin costs for simplicity and any portfolio effects. 

D. Pricing with Stochastic Severity: Non-Linear Instruments 

To price options we require stochastic intensities or severities in the Dirac process which we 
now consider. There are many possible models and we illustrate a simple and direct approach for 
stochastic severity, the Dirac-OU-Severity model. 

In all models the key is information flow and in this regard a Dirac process poses a particular 
challenge because it is uniquely memoryless. When there is no Dirac spike the process is zero. 
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Severity 



Figure 9. Illustration of Equation 31 showing b = 


where the driver maps to [3/4,1], see text for details. 


1 where the driver maps to [0,1] and 5 = 3 


This is unlike diffusion or jump processes which store some state in their current value. In a Dirac 
process this information can be stored in many places but ultimately appears in either the spike 
intensity, or in the spike severity or both. This stochastic severity or intensity is driven by an 
external process, just like the jump size for a jump process. 


Example model: Dirac-OU-Severity The Dirac-OU-Severity model uses a tanh-transformed 
Ornstein-Uhlenbeck process for the severity with a deterministic arrival intensity for Dirac spikes: 


Severity driver dx = 9{fi — x)dt + adB (30) 

Severity v = ^ + (5 _ i )/5 ( 31 ) 

2x 

= A+("-')/'> (32) 


Equation means that the transition density function is known analytically. It is conditionally 
Normal with variance and mean as below 


Var(OU)(0, fj, t) 


^2 (1 _ 
20 


Mean(OU)(a; 0 ) =h (^1 “ 6 


(33) 

(34) 


Equation 31 is the severity driver transformation, this transforms [— 00 , 00 ] 1 —)• [1 — 1/5,1], the 
severity of the event as shown in Figure]^ One minus the severity gives the probability of survival 
given an event. 

We operate in event time for the severity process, i.e. we take a scaling of unit event time 
between events. This is a common way of thinking about discrete event processes, like jump 
process, in financial modelling. Subordinators generalize the idea of event time, and this model is a 
simple example in that all times between events have the same calendar distribution (Exponential). 

In this Dirac-OU-Severity model the severity is independent of the intensity because the intensity 
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is constant. 

With this modelling setup we can price efficiently by discretizing the state space of the severity 
process. 


D.l. Zero Coupon Bonds 

THEOREM 8 (Bonds: Dirac-OU-Severity Hazard with Continuous Short Rate): The price at time 
t of a defaultable zero coupon bond maturing at time T with unit notional, assuming 

• the short rate is continuous; 

• the hazard rate is given by 


dX{u) = s{u) dD(^;W){'>j) 

(35) 

s{u) = (tanh(x(n)) + l)/2 

(36) 

dx{u) = 9{fi — x{u))du + adB 

(37) 


where ({u) is deterministic; 
is given by: 


P{t, T) = P{t, T) h (I - S)rP{N{j^ C{z)dz)) = I 


and the price of a semi-defaultable zero coupon bond is given by: 


P{t,T,,T 2 ) = P{t,n) h (^^{A{I-S)yP{N{j\{z)dz) = i}^ I 


where: A is the transition matrix; S is the diagonal severity matrix; h is a row vector expressing 
the state at t (“here”); 1 is a unit column vector. 


Proof. Proof: We operate in event time, i.e. there is unit event time between events. Therefore the 
transition matrix A is constant for each new event. At each event there is a probability of default, 
which we can put on the diagonal of the severity matrix S. By assumption the severity is driven 
by an Ornstein-Uhlenbeck process which is Markov so S is constant. The combined probability of 
moving to a new state of severity and survival of the event is A(I — S) and this is also constant. 

Since the short rate is continuous the riskless bond P{t, T) is independent of the defaultable 
part. Now in the interval (t, T] the number of Dirac events is independent of the intensity and so 
is given by a Poisson distribution with scale C{z)dz, hence 


P{t,T) = P{t,T) h Ei ^{A 

= ^ C{z)dz)) = I 


(38) 

(39) 
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The semi-default able bond price follows directly. 


□ 


Note that in general we have complete freedom to chose A providing it is a valid transition 
density matrix. Using an Ornstein-Uhlenbeck process is simply a method of generating A with a 
small number of parameters. 

We now turn to modelling the short rate with a similar process to the hazard rate. In this case 
the severity need not be limited to [0,1] but may have an arbitrary range. For generality we would 
introduce a scaling factor into the range, but for simplicity of exposition we do not include it. 

THEOREM 9 (Bonds: Dirac-OU-Severity Hazard and Mixed Short Rate): The price at time t of 
a defaultable zero coupon bond maturing at time T with unit notional, assuming 

• the Dirac part of the short rate and hazard rates are given by: 


dr{u) = sidDc_^{Wi{u) + SadD(^^{^){u) 

(40) 

d\{u) = S 2 dF>^ 20 )(^) + 

(41) 

Ci Cj: * 7^ J) 

(42) 

and 

(43) 

P{t,T) = Pc{t,T)PD{t,T) 

(44) 

where Pc{t, T) is the part of the price given by a continuous short rate and Poit, T) the part 

of the price given by a Dirac process 

the severity of the three Dirac process are given as below and use 
transformation, 

Equation 31 for range 

dXa = OaiHa “ Xa)dt aadBo 

(45) 

dxb = Obifib - Xb)dt -1- abdBo 

(46) 

dxi = — xi)dt + aidBi 

(47) 

dx2 = 02{p.2 — X2)dt + a2dB2 

(48) 

dBidBj = 6ijdt, i,j £ {1,2,3} 

(49) 

Xi{t = 0) = 0 i £ {a, b, 1, 2} 

(50) 


then the price of a semi-defaultable zero coupon bond is given by: 


P{t,Ti,T2)=Pc{t,T)llhk (f2{AkiI-k)yP{N{r\kiz)dz)=i}] I (51) 

kea \i =0 J 

a = {a,b,l,2} (52) 

Tq = {Ti, r 2 , Ti, r 2 } (Used in sequence with a). (53) 


where: A^ is a transition matrix; is a diagonal severity (credit) or change (rates) matrix; hk are 
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a row vectors expressing the state at t;l is a unit column vector. 


Proof. Proof: The non-Dirac part of the riskless bonds has no jumps so is independent of all 
the Dirac terms. This is now just a generalization of Theorem so we only need to show that the 
elements in Equationare independent. Note firstly that the driving processes for the severity are 
independent by construction, Equation 49 Also the Dirac process in Equation 42 dD(^., i G 1, 2, 3 
are independent by construction. □ 


Since we have the price of defaultable zero coupon bonds we have the price of defaultable coupon 
bonds. Since we have the price of semi-defaultable zero coupon bonds we have the price of CDS 
from Theorem m 


D.2. Option Pricing 

Pricing options using Dirac processes is simple, in the Dirac-OU-severity model using the setup 
above, because we have discretized the state space. It is further simplified because the different 
driving processes are independent although they may contribute to both credit and rates parts of 
defaultable bonds. 

Considering Theorem and Equation pricing an option on a bond is obvious, but long, 
because we already have the vector u giving the current state in the equation. We also have the 
future payoff of I (one) in all states of the world. Thus, since we have state-dependent semi- 
defaultable bond prices we automatically have state-dependent CDS prices and hence options on 
CDS. 

THEOREM 10 (Bond Option: Dirac-OU-Severity Short Rate): Let the short rate be: 


dr{t) = s{t)dDi^(^t){t) (54) 

s{t) = (tanh(x(t)) -|- l)/2 (55) 

dx{t) = 9{fj, — x{t))dt + adB (56) 


then a European call option with exercise date Tx and strike K on a riskless zero coupon bond with 
maturity T is given by: 


fTK—t 


CallZCB{t, K, Tk,T) =h ^(A (I - S)yP{N{ / Ck{z)dz) = i] 


,i=0 


rT-Tii 


max 


J2{A{I-S)yP{N{ Ck{z)dz) = i}l-K,0 


\i=0 

t <Tk < T 
K>0 


(57) 

(58) 

(59) 


Proof. Proof: Both A and S matrices are constant by construction. The max in the second line of 
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Equation |57| above gives a column vector of payoffs depending on the state of the world. Then the 
first line discounts this back to t, and the h row vector obtains the price in the current state of the 
world. □ 


Theorem 10 shows that option pricing with a single factor has the same complexity as bond 
pricing. In contrast when more factors are involved the option payout links the factors on the 
exercise date and produces a higher complexity. 


THEOREM 11 (Bond Option: Dirac-OU-Severity Hazard and Short Rate): With conditions as 
for Theorem^ except that we take Pc{t,T) = 0, the price of a European call option on a semi- 
defaultable zero coupon bond with strike K, exercise date Tk is 


CallZCB{t,K,TK,TuT2) = h | ^^^pi{TK,i)p2{TK,j)paiTK,k)pbiTK,k) 

i j k 

x'EEE max(pi(ri -TK,i)p2iT2 -TK,j)pa{Ti -TK,k)pb{T2 -Tx^k) - K, 0)1 | (60) 

i j k 

PaiT,i) = (H„ (I - 5„))*P{iV( r CfEa){z)dz) = i] 


a G {1, 2, a, b} 
/3(a) = {1,2, 0,0} 


(61) 

(62) 

(63) 


Proof. This is a direct generalization of Theorem 10 


□ 


□ 


Given Theoremj^and Theorem 11 it is clear that we can price options on CDS with the different 
levels of complexity depending on which processes we make deterministic and how many driving 


processes we use. We can assume as in (Brigo and El-Bachir, 2010) (based on Brigo and Mercurio 


(2006)) that the volatility of riskless rates has little effect and model these as constant. This is 


supported here under the assumption that the short rate is continuous. In this case there is zero 
correlation with the Dirac process for the hazard rate. On the other hand if we assume that rates 
are Dirac and that there are significant common shocks then the assumption of independence is 
false by construction and there may be significant interaction. 

In any market we may see significant common jumps between large company defaults and 
interest rates where the company is perceived to have a structural significance. Equally, if the 
government shocks the interest rates this may cause a shock to many company integrated hazard, 
i.e. survival. Thus it is a matter of calibration, and judgement, how to model the situation 
in general. In the interest of space we will only consider the simplest case to demonstrate the 
capability of Dirac processes. 

Few previous models have been able to deliver the very high implied volatilities, around 100% 


required in the CDS option space (Stamm et al., 2015). Here we will show that Dirac-based models 


can deliver this level of implied volatility. By implied volatility we mean the volatility backed out 


23 













from prices using CDS option market models (Brigo and Mercurio, 2006). These models are direct 
analogies of Swap Market Models in the interest rate space. 


THEOREM 12 (CDS Option: Dirac-OU-Severity Hazard Rate): Let the short rate be deterministic 
and the hazard rate as before a European call option with exercise date Tk and rate R on a CDSTa,Tb 
with maturity T is given by: 


CallCDSafi = h E [ P{0,TK){ProtectionLega f,{TK) — PremiumLegg^i,{TK,R))~^] 


where the protection and premium legs use the future value of the semi-defaultable bond 

p(ri,ri,r2) = {i- s)yp{N{j\{z)dz) = i}^ i 


As before underlines indicate vectors, and these are column vectors above. Underlined operators, 
i.e. max , operate element-wise. 


Proof. Proof: Now P{Ti,Ti,T 2 ) is a vector where the entries correspond to future states of the 
world. The result thus follows directly from the definition of the semi-defaultable bond and the 
discretization of the state space. □ 


Note that pricing an option on a CDS requires marginal effort beyond pricing an option on a 
zero coupon bond. This is a direct result of our state-space discretization approach. 


D.3. Example CDS Swoption Pricing 


We will now consider an example of CDS swaption pricing to illustrate the high implied volatility 
possible. We first note informally that CDS spreads generally move mostly in parallel (as for interest 
rate yield curves) except when reference entities are distressed when their CDS spread curves invert. 
When CDS spread curves invert liquidity moves from the 5Y CDS contract which is usually the 
most liquid to the lY CDS contract and data on longer-dated contracts should be viewed carefully 
to ascertain whether it is reliable (i.e. executable in comparable volume). This observation on 
parallel CDS curve moves suggests that there may be a significant term structure of volatility 


(TSOV) as in Figure 10 We use a piecewise-flat TSOV with two sections, up to option exercise 


date and beyond that date (ui, < 72 ), to obtain CDS swaptions with the implied volatilities shown 


in Figure 11 Example CDS swaption parameter settings are below: 


• 6 = 6 , tstep = 1) = 0.73, fJi = 60%, (72 = 10%, 9 = 0.1%, recovery = 40%, r = 2%, 

intensity= 2 . 

• r is the float zero yield curve level. 

• tstep states that unit event time occurs between events for the severity-driving OU process, 
and the intensity means that two events are expected to occur per unit calendar time. 
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Figure 10. Term structure of volatility shape suggested by observation of mostly parallel moves 
of CDS curves, i.e. 1/y/T. 


CDS option, 1 into 5 
CDS ATM =200bps 



Figure 11. CDS swaption implied volatility smile for l-into-5 CDS swaption example. 


This model is quite parsimonious with a relatively small number of parameters (six) for the Dirac 
part. 


III. Discussion and Conclusions 

In this paper we have introduced Dirac processes, based on sequences of Dirac delta functions, 
to mathematical finance, focussing on short rate (interest rates) and hazard rate (credit) models. 
Dirac processes provide short-rate-type models with the same expressivity that forward-rate (either 
instantaneous, HJM, or tenor, LMM) already posses when these include jumps. In fact every for¬ 
ward rate model with jumps implicitly defines a short-rate-type model with Dirac delta functions 
so these models have always existed. Thus it is now possible to choose the most appropriate model 
setting, short-rate, instantaneous forward rate, or tenor forward rate, without compromising prac- 
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tical expressivitjj^ In fact from an engineering point of view short rate models mixing continuous 
and Dirac processes can be expected to be sufficient as jumps are qualitatively redundant. 

We developed pricing using Dirac processes for basic interest rate and single-name credit prod¬ 
ucts, both linear (zero coupon bonds and swaps) and non-linear (options). Dirac processes applied 
to pricing CDS swaptions demonstrated that this approach could produce high implied volatilities 


problematic with other approaches ( 

Jamshidian 

2004 

Brigo and Mercurio 

2006 

Kokholm and 

Nicolato 

2010 

Brigo and El-Bachir 

2010 

Roti 

2013 

Weekend 

2014; 

Stamm et al. 

2015 

). In 


general a process Pt with diffusions Mt, jumps At, and Dirac processes Dt can be written as: 

Pt = Mt + At + Dt. 

However, in short rate modelling the fact that jumps add no qualitatively different behaviour (after 
integration) suggests that jumps may be redundant whilst Dirac processes expand the expressivity 
of the short-rate, and hazard-rate , modelling paradigms. 

Requiring Dirac processes in a short rate setup means that there are jumps in forward rates 
and also in swap rates. Thus Dirac process pricing potentially provides an additional tool for the 
detection of incomplete markets. However, this approach only requires a single implied volatility 
unlike methods which require a whole volatility smile to fit (i.e. using tail shape information). Given 
that the historical (and implied) CDS volatilities are high enough to make Dirac process modelling 
attractive this suggests that there are jumps in CDS swap rates and hence that the market is 


incomplete. Modelling using integrated survival probabilities (Peng and Kou, 2008) also required 


jumps, supporting this implication. Now if a market is incomplete then no option hedger can be 
profit-and-loss (PnL) flat in all states of the world (i.e. all outcomes). Any rational participant will 
require compensation for potential losses, and PnL volatility. This moves pricing into a combined 


real-world and risk-neutral pricing setup (Kenyon and Green, 2015) or, equivalently, specialized 


measures (Cont and Tankov, 2003). The required compensation may be higher than the market is 


willing to pay, thus restricting liquidity as observed Carver (2013). Of course liquidity can be low 
for many other reasons as well. 

Dirac processes are a subset of Generalized Processes that have not previously been directly used 


in derivative pricing (Seydel, 2012) because the dollar value of a Generalized value is meaningless. 


For example what is the dollar value of a Dirac delta function at the origin? This is not even a value 
in M. Generalized processes have a natural connection with stochastic partial differential equations, 
but this is a separate research area. Stochastic backward differential equations are different again. 

Any derivative modelling framework can use Dirac processes, either directly in short rate models, 
or via a time-integral approach otherwise. The method to include them depends on whether the 
underlying being modelled is a tradable or not. Short-rate models are based on a non-tradable 
that must be integrated (at least) before it becomes a tradable so Dirac processes fit naturally. 
For models of tradables the Dirac process must be integrated, possibly over very short periods. 


^Whilst HJM models are theoretically infinite dimensional no practical implementation has more than a very small 
number of factors. 
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before incorporation. For example Dirac processes fit naturally into power pricing for modelling 

short interval contracts which can go down to five-minute slots. Non-storable commodities, or 

commodities with no local storage naturally exhibit spike behaviour. 
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